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Abstract 

We compare several definitions of tiie density of a self-bound system, 
such as a nucleus, in relation with its center-of-mass zero-point motion. 
A trivial deconvolution relates the internal density to the density defined 
in the laboratory frame. This result is useful for the practical definition 
of density functionals. 

Recently [1], Engel pointed out that the usual form of nuclear Hamiltonians, 

i=l i>j = l 

forbids a proper definition of the nuclear density. Indeed, the center of mass 
(CM) of the system delocalizes in a plane wave; the nucleus is everywhere, the 
density is flat. The conclusion holds whether the two-body interaction v is local 
or non local and still holds if three-body forces are also present in H, as long as 
any explicit density dependence is avoided; Galilean invariance must be ensured. 
An elementary modification of H, 

n = H + Auj^ R'^/2, R = A-^J2^i^ (2) 

i 

traps the CM. The ground state of H now factorizes as a product of a Gaussian 
for this CM and an "internal" wave function of the {A — 1) Jacobi coordinates, 

5'(ri,r2, ...,rA) = T(R) ■ipint{^i,£,2, ••■,Ca-i) , 

. ^ r2+ri rA-i+rA-2 + ■■■ + ri , , 
£,i^r2- ri, ^2=^-3 ^ , ••■ =rA A~l ' ^ ' 

Calculations in the {ri} representation are much more convenient than those in 
the {i?, ^j} one, for obvious symmetrization reasons. The laboratory density, 

p{r) = A J dridr2 ...drA-i |*(ri, r2, r^^i, r)p , (4) 
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is much easier to calculate than the "internal" density, 

(7(0= A J d^,d^2.-dU-2 |V^„t(a,6,-,?A-2,0l', (5) 

Throughout this note, wc shall use the word "internal" instead of "intrinsic" 
when we refer to properties independent from the CM. This is because we retain 
the word "intrinsic" for those states out of which rotation bands and/or parity 
vibrations are modelized. Our understanding is that the adjectives "internal" 
and "intrinsic" belong to completely distinct concepts and models. 

Three remarks are in order at this stage, 
i) The density of interest for a density functional theory (DFT) is p, not a. 
Indeed the Hohenberg-Kohn theorem derives from embedding the system in an 
external field, namely replacing H by 

A 

lC = n + Y,u{ri), (6) 

i=l 

then considering the density T(r) of the ground state E! of /C, 

T{r) = A j dridr2 ...drA-i |H(ri,r2, ...,rA-i,r)|^ , (7) 

and finally proving that there is a one-to-one map between u and t. The ground 
state energy of K. receives the contribution / dr u{r) T{r), out of which u and r 
are recognized as conjugate Legcndre coordinates, hence the functional Legendre 
transform which defines the density functional. When u vanishes, p is that limit 
of T which minimizes the functional. It seems obviously very diSicult to set any 
similar chain of arguments in the Jacobi representation! 

ii) Despite this priority of p for a DFT, the internal nature of a is com- 
pelling for a physical interpretation. It is tempting to calculate a in the {rj} 
representation. With the Gaussian F square normalized to unity, one finds 



/ 
/ 



A-V(0= J dRd^id^2.:dU-2 [F(i?)]2|V',„t(Ci,6,-,?A-2,e)l' = 
dRd^t...d^A-2d^A-l [nR)? \i>int{^U^2,...,U-2,U-l)? = 

dn...drA-idrA S (rA - !l±lL±£_±Iizi _ |*(ri, r2, r^)p. (8) 



iii) the last Jacobi coordinate also reads 

c -4rA - (ri +r2 + ... + ta-i + va) A 

The density a therefore also represents, except for a trivial rescaling factor, the 
density referring to the internal degree of freedom, va — R- 
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Finally, according to Eqs. ^ and the density p reads, 



^ V('') = dridr2 ...drA-idrA 6{rA - r)\^{ri,r2, ■■■,rA-i,rA)\'' 



J dRdi,...dU-is(^R+^-T^U-i~r] [r(i?)]2|^,„,(^i,...,^^_2,^^_i)l2^ 



A 



J dRd^,...dU-2 [r(i?)]2 

A-^ J dR [T{R)f a 



£.1, ■■■,£,A-2, 



A 



A-1 



(r-R) 



A 



A-1 



(r-R) 



(10) 



The convolution transforming a into p is transparent, with again an inessential 
rescaling factor A/ (A — 1). The zero-point motion of the CM blurrs the internal 
density in a way which can be easily inverted, via a deconvolution. The bottom 
line is, p contains the same information as a. It is thus possible, and likely much 
easier, to design a DFT with laboratory densities t. At the end one recovers the 
internal ct by a deconvolution of p, that solution obtained by the minimization, 
with respect to t, of such a "laboratory Hohenberg-Kohn functional" . 

It will be noticed that the presence in /C of external potentials u{ri) cou- 
ples the CM degree of freedom R and the internal ones There is, in general, 
no CM factorization for eigenstates of /C. The factorization occurs at the limit 
u = 0. Then one must verify that the Fourier transform of p shows the factor- 
ized, Gaussian decay at large momenta, implied by the convolution, Eq. ()10|) . 
Otherwise, deconvolution will fail. This "deconvolution syndrome" is very well 
documented in the literature about generator coordinates. For the DFT, expan- 
sions of the density in harmonic oscillator functions and related polynomials, 
constrained [5] to satisfy matter conservation, make a useful precaution to avoid 
the deconvolution syndrome. 

For atoms and molecules, CM traps factorizing other wave packets F than 
Gaussians might be convenient, but the link between a and p remains the same. 

We conclude by claiming that a density functional theory for self-bound 
systems is available with densities in the laboratory system, without any loss of 
information about the internal, physical density. 

Acknowledgement: The author is indebted to B.K. Jennings for calling his at- 
tention to the work of Engel . 
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